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XIX. On a New Method of Approximation applicable to Elliptic and Ultra-elliptic 
Functions. — Second Memoir*. _By Charles W. Meeeifield. Communicated ly 
W. Spottiswoode, ^sg'., F.B.S. 

Eeceived March 20,— Read April 3, 1862. 

Since my first memoir on this subject was read before the Society, Mr. Sylvester has 
published a method, more general than mine, of applying rational approximation to 
facilitate the computation of the integrals of ^rational functions. This method, at 
which he had arrived independently, included, a majori, the one which was the subject 
of my memoir. Aided by his papers, my subsequent studies have enabled me to view 
the method with more generality, as well as with more precision and completeness of 
detail, and I am now able to present it in a sufficiently finished and practical form for 
the immediate use of the computer. I have also computed auxiliary Tables, to render 
its application easier in certain cases. 

Any rational formula, which gives approximately the value of a function to be inte- 
grated, may be integrated in lieu of it, and the result will in general be an approximate 
value of the integral sought. But for such a process to be of any practical utility, the 
convergence of the formula must be excessive, for the complexity of the integral forms is 
so great that the labour would be enormous, unless the terms were very few in number. 
In the discovery of formulae sufficiently convergent for the purpose, lies the success of 
the method. 

We are by no means restricted to functions under a square root, or even to pure 
radical forms at all. The principle applies with equal generality to functions which are 
given implicitly as roots of equations, and thus to a class of differential equations ; and 
Mr. Sylvester has well remarked that these formulae not only afford facilities for com- 
putation, as by a method of quadratures, but also enable us to assign superior and infe- 
rior limits to an integral, without losing its generality of form. 

I shall begin with the approximation to the square root, giving it in its general form, 
and explaining its exact analytical signification. I shall then show its application to 
Elliptic Functions, and how, in the ordinary cases, certain simple reductions can be 
effected, which greatly lessen the labour of computation ; and I shall give these reduc- 
tions for the cases more commonly occurring, with some examples and working formulse. 
I shall then add a short account of the extension of the method. 

The paragraphs in the first two sections of this paper bear a consecutive number foi' 
convenience of reference. 

* For the First Memoir, see the Philosophical Transactions for 1860, p. 223. 
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Section I. — Approximants to the Square Soot. 

1. Mr. Sylvester gives, for the approximants to the square root, the following state- 
ment : — 

" Let r be an approximate value of v^N ; then by that mode of application of Newton's 
method of approximation to the equation ^'■'=N, which is equivalent to the use of con- 
tinued fractions, we may easily establish the following theorem, viz., that 
jf + N r^ + SrN r* + 6r^N + N^ r^+I0r^N + 5rN^ 

will be successive approximations to \/N," 

2. Their general form is 

2'~(r+VN)*-(r-v'N)*^^^' ^^O 

which is always rational. In this form the approximation to v'N as i increases is 
obvious. The method of my previous memoir is simply the particular case of i=2*. 

3. If we wish to approximate to N~*, we may take the reciprocal of (I.), or, what is 
simpler, we may divide (1.) by N, thus obtaining 

(r+^N)'+(r--/N)' 1 .^ . 

Before we can integrate these formulse, we must reduce them by means of the method 
of rational fractions ; the simplest and most general way is as follows : — 

4. Let § be an ith root of unity ; then, obviously, 

log (l-x')=log{l-gx)-]-\og{l~fx)+ +log(l_/a^). 

Multiplying the differential coefficient of this by (—x), we obtain 

ia?' gx ^x ^x , f^ ■ 

1 —x*~ 1 — gi^""" f_ g2^+ 1 _ g3^+ "TT—g'x' 

and smce j^_^,=l + ^--. and ^-^,=1+^-^,, 

l—x' l — gx'l — g^x~l — g^x~ ' 1— §V 

.l+«*_l + g« 1+g^ 1+j^ 1 + g'a? 

Making x=~_^j^, we may thus divide N='* into i fractions, each of the form 

1 (r+ ^N) +§'-(>•- -/N ), 

k being any integer not exceeding i. 

5. If we add the pairs k and i—k, we obtain for the sum of the pair, 

2(r+-/N)^-2(r-A/N)^ 

8rN 8r 

— 2 (r^ + N) - (g* + g'-*^) (r^ - N) ^^ i^E^erwinmS'or' 
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according to whether the upper or lower sign be taken. Now, because ^ is an ith. root 



2k 



of unity, §''-{-§'-'' =2 cos— ^, and the sum of the pair reduces itself, for -v/N, to 



yk-- 



4rN 



2rN 



(^2 + N) - cos ^ (r^- N) * r^ sin^ ^ + N cos^ ^ 



2rN 



, kir 



2rN 



kir 



N + sin2— (»-2-N) r^-cos^-^- (r-^-N) 



(30 



For N * we have the simpler forms, 



Zi.=r 



4r 



* (^2 + N)-co8?^(r2-N) 



2r 



,kii 



,kic 



2r 



r2sin2-T- + Ncos2^ 
2r 



N + sin^— (r^-N) * r-^-cos^^ (?-^-N) 



(4.) 



All that remains is to integrate these terms, and sum them. 

6. Our gi'ouping the terms in pairs has limited the value of Jc to range from ito^{i—l) 
when i is odd. There is an odd term which, however, presents no difficulty, being 

simply - in the case of -/N, and -^jv in the case of N~*. When i is even, k is limited to 

r® + N r^ + N 

range from 1 to ^i — 1, and the odd term becomes — -■ — in the case of \/N, and —^^k— 



in the case of N~*. It is important to bear in mind that the term just mentioned is an 
odd term, and therefore not affected with the coefficient 2, which appears in the terms 
composed of pairs corresponding to imaginary roots. 

7. The value of i, which I consider to be most useful for general purposes, is i=8 : in 

this case the odd term becomes — g~- or g->j ■, and the other values of -j are three in 

number, viz. 22° 30', 45°, 67° 30'. With proper precautions i=8 will almost always 
give seven or more figures correct. 

8. If we now give infinite values to k and i and pass from the summation to the defi- 



nite integral, we have (putting A=t j 

Jo r* + (N-?-2)cos2Aw' 



:N-i= 



2rd\ 



2 r C^" d(p 

'^'> 1+ „ sm^4> 



and since 



Jo ^ 



d<p 



h 



this is an identical equation, as it ought to be. 

9. This use of approximants, therefore, is simply the application of the method of 

3l2 
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quadratures to a definite integral, which we substitute for the surd proposed for 
evaluation. 

10. It would appear at first sight that a full application of the method of quadratures 
in the ordinary way, with the help of difierences, would give better results than the 
mere summation of the ordinates. But this is not the case ; for the differences diverge 
immediately. If we use difierential coefllicients for the quadrature, instead of difierences, 
we have an opposite anomaly, namely that the correction of the summation appears to 
be absolutely nil, inasmuch as the difierential coefficients which appear in the series are 
all of odd order, and the numerator of each of them contains the factor sin (p cos (p, which 
vanishes at both the limits and ^t. Legendse has discussed this point. See the 
Appendix to the second volume of his ' Fonctions Elliptiques,' p. 578. 

11. The application of the method to integrations, then, lies in the substitution for 



I 



'« M ,, „ rT* 2Mr.dk. dt 
at 01 



in which, since X and t are perfectly independent of each other, we may change the 
order of integration, thus obtaining 

2Mr . dt 



Jo I Jo'- 



d\ ; 



and the rest of the operation depends upon our being able to perform the integration 
in } generally, and then to determine the integral in X by quadratures. The great 
advantage of the method turns upon the easy application of the method of quadratures, 
in consequence of our not requiring to difference the ordinates. 

12. One way of exhibiting generally the degree of convergence is as follows : 
N ** always lies between 



(r+ ^/N)«— (r- v'N).- '^'^^ ^^ (rH- 'v/N)'+ [r- -/N)*' 

and the error of either is therefore always less than their difference, 

^M 4(r^-N)^ 

^^ (r+v'N)^*-(r-'/N)2*" 

13. There is another mode, by which, in any given case, we may see how far it is 
necessary to carry our work in order to obtain a given number of decimals correctly in 
the result. Let 6^ be determined by the equation 

/N r 

and let sin 6^ — ^^-;^ or — tj^, whichever may be less than unity ; then the rath approximant 

/N 
will be ^~«"- This is easily seen from the general term of the approximant, since 






APPLICABLE TO ELLIPTIC AND ULTEA-ELLIPTIC FUNCTIOlSrS. 



421 



14. A table of meridional parts, such as is given in the books on Navigation, if carried 
far enough, vi^ould solve this equation. I have calculated an auxiliary Table for the 
purpose, as foUow^s : — 



Let cosec^— 1=2, log, tan (-T-{-x(p\=9/; then 

2'=ilog, ^=|log2-^logz+ilog (l+|z) 
=^log2-^log2+|-fg+i8-^+ 



To bring this formula to the same unit as the common Table of meridional parts, 
we must multiply it by the number of minutes in the arc equal to unity, or by 
L=3437-74677 07849 4, whence we have i L log, 2 = 1191-43224 08243 2, and 
iLlog, 10=3958-85223 39129 100. These data give the following Table, the argument 
being the common logarithm of z with its sign changed ; theit is, the number of places 
which are correct: 



-logz. 


y- 


1 

2 
3 
4 
5 


6234-14859 

9117-70966 

13068-84816 

17026-92712 

20985-70200 



-log«. 


y- 


6 

7 

8 

9 

10 


24944-54650 
28903-39796 
32862-25012 
36821-10235 
40779-96458 



-log z. 


y- 


11 

12 
13 
14 
15 


44738-80681 
48697-65905 
62666-51128 
56615-36362 
60674-21575 



15. As a simple example, let N=3, r=2 ; 



V 



/N 



= sin 60°: the meridional parts for 



60°= 4527 ; and in order that the error may not exceed unity in the tenth place of 

40780 
figures, we must have m or i=-r^Y = 9 ; so that we must make i=9 at least, for the 

10th figure to be correct. 

16. These methods of course only exhibit the degree of approximation on the surd 
itself. The proportionate approximation is generally greater on the integral than on the 
simple surd, because the first approximant is usually so chosen as to be identical with 
the surd at one of the limits, and it is only near the other limit that the discrepancy 
tells. 

Section II. — Details of Seduction and Computation. 

17. The chief assistance, which can be provided a jpn'on for the computer, consists in 

r* 2M.rdt 

the exhibition and discussion, for the ordinary forms, of the integral \ j^ m_j.^\ ~^;fy^j;^ 

and of the auxiliary functions which present themselves in its reduction. 

18. In applying these methods to elliptic integrals, the radical and the first approxi- 
mant r must both be of a simple form, and it is advisable thatr^— N or "N—r^ should be 
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of a square form. For the common form of the elliptic radical v'(l— sin^ ^ . sin^ ip), our 
choice is practically limited to 

(1) r=l, (2) r=cos<p, (3) r=siii ^ . cos ^, 

(4) r=cos &, (5) r=cos 6 . sin <p. 

And on these suppositions I now proceed to the integration of the general form of the 
reduced approximant for ^^(l—sin^ d.sirv' ip)~^d(p:=^zd<p. I omit mention of the con- 
stants of integration, because very slight changes in the function may alter them. The 
first of our three cases require, as they stand, no constant, and these are the most useful 

cases. 

(1) r=l, r^-N=sin'^.sin>, 

_____ 2 

1 — sin^ 9 . cos^ -.- . sin^ a 



Jzj(??)=2/^l-sin=^ . cos'y^ *tan-'|(^l-sin^ 6 . cos' yVtan®!- 

(2) r=cos(p, r*— N=:— cos''^. sin'®, 

2 coS(p 



2 cos if 



i 



cos'"* ip + cos^ fl . cos^ ^ . sin^ <p 1 — ( 1 — cos® fl . cos® -;. ) . sin® ^ 



t / 



^^dip = ( 1 — cos' ^ . cos' — j log, 



l+sinipi 1— cos®d . cos® -^1 



k7t\i 



1— sinip(l— cos® 9. cos® -A 



(,3) r=sin 6 . cos <p, r'— N= —cos' 0, 
2 sin 9 . cos <p 



Zu — 



A - cos® 9 . sin® ^\ - sin® 9 . sin® ^ 



^z^dip = ( 1 — cos' ^ . sin' -j- j log - 



{ 1 —cos® 9 . sin® -r- 1 +sin 9 . sin <p 
{ 1 — cos® 9 . sin^ -^\ —sin 9 . sin f 



(4) r=cos^, r'— N=— sin'^cos'(p, 

2 1 

1 + tan® 9 . cos® -r- . cos® ip 



^Zt:d<p=2 ( 1 — sin' 6 . sin' j\ ' tan-'j cos . tan ^/ 1 — sin' 6 . sin' j\ ' I. 



(5) r=cos^sin®, r'— N~— cos'?), 

2 cos 9 . sin f 

cos® 9 . sin® ip + cos® -;- . cos® f 
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j-M9=2 (cos^ ^-cos^ ^) '* tan- j^ (cos^ ^-cos^ ^) "*} 

cos 9\ I ) 



= 1 cos^^ — cos^^-^ j logj 



H ^(cos*S — cos^ -7- 1 

cos 9\ * / y 



(6) If we make i=tan ^^, we obtain 

(1-sin^ 6 . sin' ^)-*#=2(l-2 cos 2 ^ . f +#')-* cZ;^. 
Taking r=■\—^^, the terms which we have to integrate are of the form 

J(l-i!2)H4cos2fl . cos^^ . <2 
Putting §'=1 — 008' 9 . cos' -j, we have 



J.,#=-Jlog.0^|±-J). 



The same expression serves for the integral 

Idt 



i 



V'(H-2cos2fl.if2 + if«)' 



if we put §-'= 1 — sin' ^ . cos' — . 

19. It win be observed that the first four cases, and the sixth, depend upon a radical 

of the form -/(l— sin' A. sin'<y), where u is restricted to the selected values of ~. 

Assuming the modulus sin A not to vary, it would therefore in general be better to begin 
by computing the radical for the selected values. I have computed, and I append to 

this paper, a Table of this radical, the selected values of y being 22° 30', 45°, and 67° 30', 

while A ranges by whole degrees from 1° to 90° inclusive. Every entry but the last in 
the 2nd, 3rd, and 4th columns of the Table was computed by myself in duplicate with 
Vega's ten-figure logarithms, by the help of two or more of the following formulae, some 
of which are from Legendre. 

20. Putting A for v'(l — sin' A . sin'a), 

(1) Make sin A . sin <i;=sin M ; then A=cos M, 

log sin M=log sin A+log sin », log A=log cos M ; or else 

(2) Make tan A . cos <y=:tan M ; then A=cos A . sec M, 

log tan M=log tan A+log cos 4*, log A=lpg cos A+ar. co. log cos M. 
Moreover, let L be the tabular angle nearest to the anglfe M: it is not necessary to 
obtain the value of M : so that we have simultaneously, 
log sin M=log sin L+s, 
log tan M=log tan L+^, 
log cos M=log cos L+c ; 
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then we shall also have, and with great approximation, 

log s=log (f . cos^ L)^(t—t . cos^ L) 

=log (c . cot^ L)+(c+c. cot^ L), 
log c=log (s . tan^ L)+(s+s • tan^ L) 

=log (t . sin^ L)±(t-t . sin^ L), 
log i5=log (s . sec^ L)ip(5— s . sec^ L) 

=log (c . cosec^ L)+(c— c cosec^ L). 

I have given the whole set of six, but my Table was computed with the pair for log c. 
By way of example, I add a specimen copy of one of my working sheets. The use of so 
many as ten figures is not altogether unnecessary, because otherwise, when A is nearly 

equal to unity, the value of log (1 — A) or of log r-— ^ cannot be had with exactness. 

21. The following formulse will also be found in many cases preferable, both for 
exactness and facility, to the ordinary use of logarithmic tables by means of differences. 
These formulae, as well as those of the previous paragraph, are but applications of 
Taylor's theorem, reduced to a shape fit for the computer. Even where ojily seven 
figures are required their application is frequently much easier, and gives more exact 
results, than interpolation by differences. In what follows, x is supposed to be thf 
nearest tabular entry. 

22. To find log ^ from log tan y. — Let us assume simultaneously 

log i/=log w+l, logtany=log tana?+l 
Putting u=logw, 2 = log tan ^, we have 

du sin 2x , d^u -.^sm2xf ^ sin 2x\ 
dz=-2^ and^,=M-^(^cos2a;-^™j, 

M being the modulus of the logarithms. 
Hence, by Tayloe's theorem, 

^=^ ~^^Y-+^*\~2^ ^^ 2^j J- nearly. 

Taking the logarithm, this becomes 

- J , /.sin 2ir\ _, /sin2a; _ \ 

log l=log[t^^-~) +t ^ -2^-cos 2wj 

, // . sin X . cos x\ t . sin x , cos x , „ . „ 

=log (^ j 4- ^ ±t+2t sm^ X. 

The latter is the better shape for a working formula, because log sin x and log cos x 
are found in the same page and line as log tan x, while log sin 2x must be looked for 
elsewhere. The first term alone is sufficient when x is small ; but when x much exceeds 
45°, cos 2x changes its sign, and even the entire formula is insufficient. The maximum 
value of the coefficient of tia the second term is 1-0631, corresponding to x=78° 33' 26"-5. 
In many cases, where the first term alone is insufficient, a rough interpolation, made at 



APPLICABLE TO ELLIPTIC AND ULTRA-ELLIPTIC FIJNCTIDKS. 



425 



sight from the following Table, will answer the purpose ; it is a Table of the value of 

(Sin ^x \ 

-~ — cos2ir I and of its logarithm, from a?=45° to ir=90°. 





46 


0-63662 


9-80387 


50 


0-73816 


9-86815 


55 


0-83149 


9-91986 


60 


0-91349 


9-96070 


65 


0-98041 


9-99141 



70 


1-02910 


0-01246 


75 


1-05658 


0-02390 


80 


1-06216 


0-02619 


85 


1-04334. 


0-01843 


90 


1- 


0- 



The Table shows that, past 45°, the formula 

, , , /^ sin a; cos x\ , 

logZ=log(^ j+^ 

is a better approximation than when the +# is omitted. It is to be remarked that t is 
at its minimum for a;=45°, and increases both towards a;=0 and a?=90°. Near the 
latter limit, where great accuracy is required, we must proceed as follows. 

Find the correction for the logarithm of the complement of the arc by the above 
process, and then find log (^ f— 3/) from log y. For this purpose, I observe that 
logy=logiF4:^ is equivalent to y=iP.10*', hence 

ix-y=ix-x.lO^'=(iT-a?)-a:(10*'-l). 

Now, let +A=10*'— 1, whence 

log (+mA)=log(4:Z)+^Z— To M? nearly, and also 



12 

+ mAa;> 



H(i,-,)=log(i.-.)-(^)-iM(|^') -. . . . 

It is not often that the third term of either formula will be required. 

I have gone into all this detail, because the inverse tangent is continually presenting 
itself in all these integrations, and because no book that I know shows the proper way 
of handling it. 

23. The following constants are needed for these and similar formulae : — 

10+logm =9-63778 43113 00537, logM=0-36221 56886 99463, 

10+logr =8-24187 73675 90828, -logl°=l-75812 26324 09172, 



10+logl' =6-46372 61172 07184, 
10+logl"=4-68557 48668 23541, 



-log 1' =3-53627 38827 92816, 
-logl"=5-31442 51331 76459. 



24. As an example of finding the inverse tangent, let it be required to find log 2^ and 
log (ir—y) from 

log tany=9-02313 50437. Here we must take 
log tan a:= 9-02303 57359 
+^=~ 
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9 93078; 



3 M 



a?=6°ri0"=21670" 
-a'= 302330" 



*7r. 
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log sin ^=9 '0206346 

logcosiP=9-9975988 

log ^=5-9969834 

ar. CO. loga^=0-9785662 

t sin 2x 



2nd correction 



5-9937830 

-14 

log ;=5-9937816 
¥ +493 



log sin^a?= 8-0413 

log 2=0-3010 

log ^=5-9970 

log(2sin^a;)=4-3393 

2sin*a?= 21,800 

t sin 2a! 



log 21670 = 4-33585 89113 

logl" =4-68557 48668 

log 0^=9-02143 37781 

1= +9 85786 

log^= 9-02153 23567 



2a? 



-= 985,8 



1007,6 

-i=-993,l 

logmA=5-9938309 2nd correction = 14,6 

log ar= 9-0214338 

ar. CO. log (|7r-jy)= 9-8339439 

log correction= 4-8492086 



The comma cuts oif the 
eighth decimal. 



log 302330=5-48048 12441 

log l"=4-68557 48668 

log (ix-x)=0-16605 61109 

correction= — 70666 



log (i^_y)= 0-16604 90443 



Verification. — ^The numbers corresponding to these logarithms of y and of ^—y are 
0-10508 29743 and 1-46571 33525, the sum of which, to the very last figure, is 
exactly -^tt. 

25. To find \o^y~^ from \o^y.— Lei log |3Y=log^4:j?, and \oq y=\.(iQ x±q; 
then log|)= log (-a^jj+rS' +"2^05 nearly. This formula obviously fails where y is 
near unity ; in this case log -^^ cannot be had with great accuracy, unless y itself be 

given absolutely. All the cases of i ^_ ^ may be included in the above formula by 

giving proper signs to p and q^. It may save trouble to remark that x must not always 
be taken to the extreme limit of the Table, because log(^+l) ^^d log (^ — 1) have also 
to be taken out. As an example, let 



log2/=0-36290 63835 
log 0^=0-36285 93030 
q= 4-4 70805 
log ff= 5-6728411 
ar.co.log(a;'-l)=9-3647540 
5-0375951 
logx=0-3628693 
log 2=0-3010300 
5-7014844 
2nd correction= —1225 



a;=2-306, a;+l=3-306, a;-l=l-306 
log(^+l)=0-51930 28492 
log(o?-l)=0-11594 31769 



sum=0-63524 60261 



-T^ =€-0001090 

iq^_ 136 

2nd corr"=0'0001225 



difierence= 0-40335 96723 
~p= -5 02761 
y+1 



log 



y- 



:0-40330 93962 



log^=5-7013619 
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This example has been so chosen as to admit of easy verification. In fact i/=2'S()626, 

w+1 529 

and log |3Y=log 2^=0-40330 93959 24. The error is therefore only of three units 

in the tenth decimal place, where there was no reason to expect accuracy. 

26. The only other formulae which I shall give are the following, for finding the 
logarithm of a number, and vice versd. They are indispensable where more than seven 
figures are required. 

Let log (x+h)=log x+k, then 

, , , /mh\ _\ mh , 

log^=log (^-j +2 -- nearly, 

log ^=log {M.xJc)+ 2 ^ nearly. 

The values of log w and logM have been given in paragraph 23. 

27. As an example of the application of the method to the evaluation of elliptic 
integrals of the third class, let us take the integral 

'^ d<^ 

^j (1— sin^a.sin^f) (1— sin^fl.sin'^^)* 

for the values a=:45°, ^=30°, (p=60°. 

I have selected these values because they can be obtained without reduction or inter- 
polation from the Table of A(^, f) which I have given, and also because sin^a=sin^, 
and therefore the integral can be reduced to one of the first class, p?ws an inverse tangent, 
thus admitting of easy verification. For this case 

2 



r 



Zic- 



■ ~ y — ~- — • 7 r"> 

(1 — sin^a . sin^ip) 1 1 - sin^ fl . cos^ -^ . sin^f | 



, , 2 sm^a 1.1/ ^ \ 

\Zudp-=- — — -;—• tan"' (cos 05. tan®) 

i 



2sin«flcos2^ 
i 



(l— sin^^.cos^YJ tan~x( 1— sin^^.cos^Y) tan?)>. 



Jcfjr 

Making Y successively 22° 30', 45°, 67° 30', and, for the odd term, 90°, we find, after a 
few obvious reductions, that eight times the value of the integral is 

{¥+A^(45°,67i)+ An45°,22i) } ^5^° tan- {cos 45°. tan 60°} 
-^^^tan-^cos 30°. tan 60°} -1 ^(30^45°) tan-'{A(30°, 45°).tan 60°} 

_ tan-' {A(30°, 67i) -tan 60°} .cos^ 22^ _ tan-' {A(30°, 22|) .tan 60°} .cos^67i 
A2(45°,67A).A(30°,67i) A«(45°, 224).A{30°,22i) 

3m2 
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As these inverse tangents range generally from 45° to 60°, I computed them by the 

(f Sin 'y oos'y\ 
'- — ) +i( ; this being sufficient to 

give eight figures of decimals accurately. I found 

Iogtan->{cos45°tan60°} =9-94747 15296, 

Iogtan-'|cos30°tan60°} =9-99246 23739, 

logtan->]A(30°,45°)tan60°} =0-00766 92607, 

log tan-|A(30°, 67i).tan 60°^ =9-99727 33807, 

log tan-' I A( 30°, 22i).tan60°} =0-01665 09657. 
I hence obtained the foUovping values : — 

For tlie positive terms. For the negative terms. 

7-10091 3039 1-13483 2441 

4-37212 6152 0-72539 4027 

2-70421 6251 1-66839 9478 



14-17725 5442 0-16728 4032 



3-69590 9978 3-69590 9978 



8)10-48134 6464 



1-31016 8183 value required 

A more exact value of the integral, otherwise obtained, is 

iF(30°,60°)+tan-'(^~|^°j = 1-31016 8161, 

which differs from the previous value by 2 units in the eighth decimal place. 

28. In order to find how many places ought to have been accurately obtained, I 

13 

observe that the method followed gives N=Tg, r=l, whence 



log /'i^N) = 9-95491 =log sin 64° 20' 30". 



The corresponding meridional parts are 5086-5, which must be multiplied by «=8, 
giving 40692-0. Eeferring to the Table in paragraph 14, I find that this nearly corre- 
sponds to ten places correct, and therefore that the integral ought to be correct to at 
least that extent. That it is not so, is due to my having curtailed the formula for 
finding the logarithms of the inverse tangents. But my object Avas only to give seven 
decimals correct, and my going beyond that was simply because, with a ten-figure 
Table, putting down the additional figures gave me less trouble (once I had to use more 
than seven) than abbreviation would have done. This remark may at first sight seem 
strange to any one who has not had some practice in using large Tables. But the loga- 
rithmic corrections are given in the shape of arithmetical complements : with reference 
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to the lOth figure, therefore, considered as an integer, the index is right as it stands, 
and we need not bestow thought on the proper placing of the correction, as we must if 
we use any other number of figures. 

29. If we had been content with five decimals, the calculation would have been very 
easy, for in that case we might have used six-figure logarithms, and have made * = 4, 
thus omitting the terms containing 22°^ and 67°^. We should get 

7-10091 3039 1-134:83 2441 

1-86022 6468 0-72539 4027 



4)5-24068 6571 1-86022 6468 

1-31017 1643 value required. 

30. It is worth while to notice a case which will sometimes occur, namely (using the 
notation of the last example), that the values may be so selected as to give, for one of 

the values of Tc, sina=sin^.cos y, and thus each of the terins into which '\^Zj,d<p was 

divided would become infinite. Of course the difiiculty is only apparent ; for in this 

case the proper value is fefe#=! (i_8inaj^sing;p)2 ' of which the integral may be at once 

,.™ ,. ,. ,, . sina.cosa 

found by differentiating the expression ^_g-^2^ g-^g. - 

Section III. — Extension of the MetJiod. 

In respect of rapid approximation and precision of limit, the foregoing processes 
leave nothing to be desired, as far as concerns the radical of the square root ; but they 
do not go beyond that. Mr. Sylvester has given an elegant extension of the method 
to radicals of a higher index, by means of symmetric functions*. 

The more general problem before us is that of approximating to the integrals of 
irrational functions by means of rational substitutions. 

Let <p and ■v^ be functional symbols, and y a function of a; then, that ^{z).y^ and 
(p{z):y^ should both be approximations to ip{z), depends upon y^ approaching unity as m 
increases. Assuming that y„ and y, are connected by the equation y„=4(m, y\), our 
problem is to choose -i^ so that, in the first place, the approximation shall be exceedingly 
rapid, and, in the next place, that <p{z).y,^ and <p{z):y,^ shall both (or at least one of 
them) be thoroughly manageable, and easily integrable. In the case of the approximants 

already given, the equation y^=-^{m,y,) has been ra&Ae ^y^—^JY::::^' 

I am acquainted with three general methods which effect the object more or less. 
The first is the obvious one afforded by the Newtonian approximation to the roots of an 
equation ; viz., let a be a first approximate solution, obtained by trial, of the equation 
/ar=0, and call f'x the differential coefficient oi fa; then a second approximation is 

* See the Philoaophical Magazine for December 1860, Supplementary Number, vol. xx. p. 525, note A. 
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^~fa~^ ; a third approximation will evidently be h—7jfj=c, and so forth. If we apply 

this method to the pure equation a;"=^, the convergent terms which we obtain are as 
follows : — 

[n-\)a''+p 

_{n— !){(«- l)a«+jo^'' + M».^.a«(»-i) 
mV-'|(w—1)«'' +;?[»-' ~' 

The second method is that of the reversion of series ; it is sufficiently discussed by 
Arbogast *. 

The third method was suggested to me by Mr. Catlet's remark that Mr. Sylvestee's 
third approximation is a particular case, for %=2, of the common form (of the books on 

the binomial theorem) v^N= /_,^^y,> ,A „ .a, approximately, « being a first approxi- 
mation. In order to gain generahty, and thereby symmetry, I shall pass from the par- 
ticular form -^/N to the more general ip~'N by the following Lemma: — 

Let N=«o+<'^i^+«2^^+«8^+ , (!•} 

and let x^, x^-, X3, x^ be determined by the system of equations, 

K=a,(l+lx) 

=«„(l+5^.(H-J^.)) \ (2.) 

=-0(1+7:^3(1+^^.(1+5..))), 
and so forth; also let N— «„=/!*, and 

l=:(l/>—aiX—a^—aiX'— )(\o+\'«+>^!^+K^+ ), • . (3.) 

then _Kq Xj a^ x„_i 

For, if we substitute these values in the equations (2.) after placing them in the 
following form, 

=aiX2+a2X2X, (4.) 

=aiX3-\-a2X3X2-\-a3X3X2Xi, J 
and so forth, we obtain 

M-^- —^-~= x^ =&c., ... ... (5.) 

which are the same equations as we should get by multiplying the two series in (3.) and 

equating to zero the coefficients of x and of its powers. The coefficient A„=-, 
obviously. 

* Calcul de? Derivations, pp. 288--296. 
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The coefficients X may now be found in a variety of ways ; by solving equations (4.) 
or (5.), by simple division, or by Abbogast's processes *. The object of the preceding 
lemma is to connect the quantities x„ with the coefficients of division and of recurring 
series. Our results in any way are, 



■ + a 



*'—«,' '^^—a^ + a,^' 



w. 



**— tt^^ + (2«g«i + Of«,2 + 3a2«iV + «/' 

If for flo) «i5 &c. we substitute the coefficients of the binomial theorem, so as to make 
N=(c+a;)", we obtain 

_ (w— l)a" + N 

__(»— l)a"+(w+l)N 
*+*2— (w+lia^ + lM-liN-^' 

(w^- l)a^" + (4w^ + 2)a"N + (»^- 1)N^ 

^+^^3— (n+l)(» + 2K»+4(n2_l)a»N^-(n-l)(«-2)N2•^• 
Making w=2, we obtain Mr. Sylvester's approximants to the square root, and A„is then 
the coefficient of «" in the development by ascending powers of 

1 

(N-^)-2aa;-a;'''' 

and so far the method agrees with the Newtonian approximation by continued fractions ; 
but from this point the two methods diverge. For w= 3, X„ is the coefficient of «" in the 
development of 



and the successive approximants are 

2a3+N «1+2N 4a6 + 19a3N+4N^ 5a9+45a«N + 30a3N^ + N3 

3a3 • ^' 2a8 + N * ^' \Q(^+ ISa^N + N^ * '^' \5a^ + 5lam + \5c^W ' ^' *^^- ' 

while the second approximant obtained by successive substitution is 

16a9 + 51a6N + 12a6N^+ 2N3 



36a9 + 36a«N + 9a8]N2 



a. 



What these methods all effect is simply a rational approximation to the value of y in 
the equation <p(y,z)=0. Then, making 5^=^, we have only to integrate in order to 
ffiid the value of u. They thus constitute a means of approximately solving, in respect 
of u, differential equations of the form ^(2, j)=0; but they do not effect the solu- 

* See his ' Calcul des DIrivations,* pp. 26, 29 ; or De Moegaw, ' Diff. Oalc' p. 331. 
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tion of this equation in respect of a, and still less do they solve the more general form 

It may suggest processes of reduction in some cases to remark, that there are many 
other functions of y^ and ^(z), which will approximate to <p{z) as m increases, besides 
the simple product or quotient of <p{z) by y^. 

There is one point about these higher approximants, of which a solution, even if 
accompanied with considerable restrictions, would be extremely desirable, — I mean the 
resolution of the denominators into factors. I do not suppose that the problem, in its 
perfectly general form, admits of a compact solution ; but any class of cases, of even 
moderate generality, for which it could be elegantly solved, would probably have very 
useful applications. The criterion of convergence and the measure of approximation 
would also have their interest. 
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Specimen sheet of work for the Table. 



Arc of 60°. 



logsin22|= 9-58283 96605 83* 
log sin A = 9-93753 06316 96 



log tan'' L= 9-0913166 
log s= 5-3160144 



log cos L=9-97473 27132 
c= 4- 25546 



log sin M = 9-52037 


02922 


79 


4-4073310 


log A= 9-97473 


62677 


log sin L = 9-62039 

s = -2 


09944 
07021 




S+-1 

V -233 
s tan* L= J 






s tan* L = 


25646 




log c= 4-4073077 






log sin 45°= 9-84948 


60021 


68* 


log tan* L= 9-7781471 


log cos L= 9-89794 


07883 


log sin A = 9-93753 


06316 


96 


logs= 4-1137763 


c= 


7797 


log sin M = 9-78701 


56338 
43344 
12995 


64 


3-8919234 

s+i 

} + 21 
8 tan* L= J 


log A=9-89794 


00086 


log sin L = 9-78701 

8 = + 






s tan'' L = 


7821 




log 0= 3-8919255 






logsin67i= 9-96561 


63459 


21* 


log tan* L= 10-2601550 


log cos L=9-77806 


24362 


log sin A = 9-93753 


06316 


96 


logs= 4-7816118 


c= - 1 


07593 


log sin M = 9-90314 


59776 
99296 
60480 


17 


5-0317668 

I +168 
s tan* L= J 


log A=9-77806 


16759 


log sin L = 990313 

s = + 






stan'L = 1 


07589 




log c= 5-0317836 






log tan22i= 9*61722 


43146 


62* 


logsin«L=: 8-6140815 


log cos g2J-=9-9666l 


53469 21* 


log cos A = 9-69897 


00043 


36 


log t= 5-5461669 


log sec L =0-00911 


84784 



log tan M = 9-31619 43190 

log tan L = 9-31615 92213 

t = + 3 60877 

tsin^L = 14018 



08 



4-1691374 



I +327 

t sinS L= J 

log c= 4-1591701 



9-97473 
c= + 



38243 
14430 



log A= 9-97473 52673 



log tan 45 =10-00000 00000 
log cos A = 



00* 



log sin^ L= 9-3009948 
log t= 6-3429218 



log cos 45 =9-84948 
log sec L =0-04845 



50021 68* 
O6OI7 



log tan M = 9-69897 00043 
log tan L = 9-69894 79790 



36 



t = +2 20253 
t sin" L = 44047 



4-6439166 
t— 
tsin^L: 

log c= 4-6439342 



::} 



+ 176 



log tan 67i= 1 0-38277 56853 
log cos A = 9-69897 00043 



log tan M =10-08174 56896 

log tan L =10-08175 16044 

t = - 69147 

tsin«L = 36075 

Ifote. — The entries marked *, 
•n tHs page corresponds to in 
MDCCCLXII. 



38* log sin2 L= 9-7730722 

36 log t= 4-7719327 

74 4-6450049 

■ I _ 24 

t sin«i L= J 

log c= 4-5450025 

and the whole of the letter-press, were printed on the sheets, 
the Table. 

3n 



9-89793 
c= + 


56039 

44049 


log A=9-89794 


00082 




log cos 671=9-58283 
log sec L =0-19621 


96605 83* 

55228 


9-77805 


61834 
85075 


log A=9-77805 


16759 



The letter A 
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Table of the value of the function log A(fl, ») or log v'(l — sin^^ . sin^ij^) for four values 

of a;, viz. 22° 30', 45°, 67° 30', and 90°. 



&. 


log A (^,22° 30'). 


log A (((,45°). 


log A (^,67° 30'). 


log 003 i. 


e. 


o 

1 


9-99999 03138 


9-99996 69274 


9-99994 35386 


9-99993 38497 


1 


2 


9-99996 12643 


9-99986 77197 


9-99977 41349 


9-99973 53589 


2 


3 


9-99991 28792 


9-99970 24074 


9-99949 17311 


9-99940 44063 


3 


4 


9-99984 62048 


9-99947 10407 


9-99909 62308 


9-99894 07898 


4 


6 


9-99975 83051 


9-99917 36910 


9-99858 74990 


9-99834 42260 


6 


6 


9-99965 22633 


9-99881 04507 


9-99796 53618 


9-99761 43489 


6 


7 


9-99952 71805 


9-99838 14325 


9-99722 96070 


9-99676 07098 


7 


8 


9-99938 31764 


9-99788 67716 


9-99637 99831 


9-99575 27764 


8 


9 


9-99922 03891 


9-99732 66254 


9-99541 62004 


9-99461 99270 


9 


10 


9-99903 89748 


9-99670 11740 


9-99433 79300 


9-99335 14689 


10 


11 


9-99883 91084 


9-99601 06211 


9-99313 48042 


9-99194 65764 


11 


12 


9-99862 09825 


9-99525 51957 


9-99183 64168 


9-99040 43940 


12 


13 


9-99838 48090 


9-99443 61615 


9-99041 23213 


9-98872 39328 


13 


14 


9-99813 08170 


9-99356 07689 


9-98887 20337 


9-98690 41185 


14 


15 


9-99785 92542 


9-99260 23560 


9-98721 50300 


9-98494 37781 


15 


16 


9-99757 03868 


9-99159 02501 


9-98544 07479 


9-98284 16370 


16 


17 


9-99726 44986 


9-99051 48183 


9-98354 85857 


9-98059 63156 


17 


18 


9-99694 18915 


9-98937 64594 


9-98153 79027 


9-97820 63255 


18 


19 


9-99660 28857 


998817 66060 


9-97940 80200 


9-97667 00664 


19 


20 


9-99624 78189 


9-98691 27246 


9-97716 82198 


9-97298 58164 


20 


21 


9-99587 70469 


9-98558 83197 


9-97478 77462 


9-97015 17377 


21 


22 


9-99549 09429 


9-98420 29331 


9-97229 68056 


9-96716 68605 


22 


23 


9-99508 98979 


9-98275 71478 


9-96968 16661 


9-96402 60827 


23 


24 


9-99467 43204 


9-98125 15899 


9-96694 41603 


9-96073 01626 


24 


25 


9-99424 46358 


9-97968 69297 


9-96408 26837 


9-95727 57115 


25 


26 


9-99380 12870 


9-97806 38852 


9-96109 61968 


9-95366 01 869 


26 


27 


9-99334 47337 


9-97638 32245 


9-95798 37258 


9-94988 08840 


27 


28 


9-99287 54524 


9-97464 57677 


9-96474 42643 


9-94593 49269 


28 


29' 


9-99239 39363 


9-97285 23906 


9-95137 67741 


9-94181 92587 


29 


30 


9-99190 06948 


9-97100 40266 


9-94788 01877 


9-93753 06317 


30 


31 


9-99139 62526 


9-96910 16705 


9-94425 34100 


9-93306 65951 


31 


32 


9-99O88 11517 


9-96714 63813 


9-94049 53211 


9-92842 04835 


32 


33 


9-99035 59484 


9-96513 92852 


9-93660 47788 


9-92359 14023 


33 


34 


9-98982 12144 


9-96308 15797 


9-93258 06231 


9-91857 42135 


34 


35 


9-98927 75363 


9-96097 46359 


9-92842 16784 


9-91336 45194 


35 


36 


9-98872 55150 


9-95881 95031 


9-92412 67607 


9-90795 76446 


36 


37 


9-98816 57654 


9-95661 79121 


9-91969 46797 


9-90234 86165 


37 


38 


9-98759 89157 


9-95437 12781 


9-91512 42488 


9-89653 21441 


38 


39 


9-98702 56072 


9-96208 12065 


9-91041 42888 


9-89060 26944 


39 


40 


9-98644 64943 


9-94974 93945 


9-90556 36388 


9-88425 39665 


40 


41 


9-98586 22425 


9-94737 76364 


9-90057 11640 


9-87777 98629 


41 


42 


9-98527 35294 


9-94496 78273 


9-89543 57674 


9-87107 34581 


42 


43 


9-98468 10429 


9-94252 19663 


9-89015 64024 


9-86412 74638 


43 


44 


9-98408 54812 


9-94004 21611 


9-88473 2O869 


9-85693 40901 


44 


46 


9-98348 75524 


9-93763 06317 


9-87916 19193 


9-84948 60022 


45 
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Table 
{continued). 



t. 


logA(<:,22°30'). 


log A {6, 45°). 


log A (((,67° 30'). 


log cos /. 


«. 


d 


9-98288 79722 


9-93498 97136 


9-87344 50980 


9-84177 12731 


46 


47 


9-98228 74657 


9-93242 18620 


9-86758 09423 


9-83378 33303 


47 


48 


9-98168 67640 


9-92982 96543 


9-86156 89166 


9-82561 08951 


48 


49 


9-98108 66053 


9-92721 57944 


9-86540 86587 


9-8I694 29168 


49 


50 


9-98048 77326 


9-92458 31150 


9-84910 00106 


9-80806 74967 


50 


51 


9-97989 08943 


9-92194 45794 


9-84264 30543 


9-79887 18039 


51 


52 


9-97929 68415 


9-91927 32846 


9-83603 81632 


9-78934 19787 


52 


63 


9-97870 63285 


9-91660 24627 


9-82928 69976 


9-77946 30249 


63 


64 


9-97812 01111 


9-91392 54820 


9-82238 76655 


9-76921 86862 


54 


65 


9-97763 89457 


9-91124 68470 


9-81634 46333 


9-75859 13013 


55 


56 


9-97696 36878 


9-90856 71996 


9-80815 89399 


9-74756 16513 


56 


57 


9-97639 47918 


9-90589 33160 


9-80083 31625 


9-73610 87646 


57 


58 


9-97583 33090 


9-90322 8IO75 


9-79337 05490 


9-72420 97077 


68 


59 


9*97527 98869 


9-90057 66154 


9-78577 61006 


9-71183 93361 


59 


60 


9-97473 52675 


9-89794 00084 


9-77805 16759 


9-69897 00043 


60 


61 


9-97420 01874 


9-89532 65788 


9-77020 61056 


9-68557 12291 


61 


62 


9-97367 63742 


9-89273 67330 


9-76224 63190 


9-67160 92909 


62 


63 


9-97316 15478 


9-89017 79885 


9-76417 74828 


9-65704 67649 


63 


64 


9-97265 94174 


9-88765 39622 


9-74601 21630 


9-64184 19615 


64 


65 


9-97216 96810 


9-88516 93618 


9-73776 04387 


9-62594 82693 


65 


66 


9-97169 30239 


9-88272 89739 


9-72943 51756 


9-60931 32999 


66 


67 


9-97123 01176 


9-88033 76506 


9-72105 11125 


9-59187 80116 


67 


68 


9-97078 16179 


9-87800 02961 


9-71262 51046 


9-57357 54170 


68 


69 


9-97034 81644 


9-87672 18497 


9-70417 63081 


9-55432 91617 


69 


70 


9-96993 03790 


9-87350 72689 


9-69572 63771 


9-53405 16846 


70 


71 


9-96952 88643 


9-87136 15099 


9-68729 96519 


9-51264 19176 


71 


72 


9-96914 42028 


9-86928 95088 


9-67892 33280 


9-48998 23640 


72 


73 


9-96877 69661 


9-86729 61579 


9-67062 76041 


9-46593 63400 


73 


74 


9-96842 76594 


9-86538 62846 


9-66244 57824 


9-44033 80750 


74 


75 


9-96809 68303 


9-86356 46269 


9-65441 43168 


9-41299 62305 


75 


76 


9-96778 49569 


9-86183 58088 


9-64657 27859 


9-38367 51767 


76 


77 


9-96749 25025 


9-86020 43167 


9-63896 37732 


9-35208 80330 


77 


78 


9-96721 99032 


9-85867 44678 


9-63163 26324 


9-31787 89102 


78 


79 


9-96696 75670 


9-85725 03968 


9-62462 71226 


9-28059 88450 


79 


80 


9-96673 58731 


9-86593 60134 


9-61799 68926 


9-23967 02300 


80 


81 


9-96652 51705 


9-85473 49940 


9-61179 28070 


9-19433 24413 


81 


82 


9-96633 57778 


9-86365 07456 


9-60606 61106 


9-14355 53039 


82 


83 


9-96616 79819 


9-85268 63874 


9-60086 74357 


9-08689 44712 


83 


84 


9-96602 20377 


9-86184 47281 


9-596"24 56795 


9-01923 46656 


84 


85 


9-96589 81672 


9-85112 82461 


9-59224 67793 


8-94029 60083 


85 


86 


9-96679 66594 


9-85053 907O8 


9-58891 24439 


8-84358 45184 


86 


87 


9-96671 73697 


9-85007 89667 


9-68627 88989 


8-71880 01636 


87 


88 


9-96566 07187 


9-84974 93212 


9-68437 67166 


8-54281 91639 


88 


89 


9-96562 66934 


9-84955 11322 


9-68322 48116 


8-24185 53184 


89 


90 


9-96561 53459 


9-84948 50022 


9-68283 96696 


— log. infin. 


90 



